Reversal Modes of Simulated Iron Nanopillars in an Obliquely Oriented Field 
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Stochastic micromagnetic simulations are employed to study switching in three-dimensional mag- 
netic nanopillars exposed to highly misaligned fields. The switching appears to proceed through two 
different decay modes, characterized by very different average lifetimes and different average values 
of the transverse magnetization components. 

PACS numbers: 75.75-(-a, 61.46-fw, 64.60.Qb 
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I. INTRODUCTION 



Nanoparticlcs with strong uniaxial anisotropy are fre- 
quently employed in technologies where the orientation 
of the particle magnetization corresponds to a binary 
digit. The switching process of the magnetization in an 
applied field therefore becomes a technologically impor- 
tant process, especially for fields that vary on a timescale 
of nanoseconds. Micromagnetic simulations have been 
essential for understanding properties of systems that 
evolve over such short time scales. 

In this paper we examine a model nanopillar where the 
anisotropy is due solely to its shape, and subject it to a 
reversal field that is highly misaligned with respect to the 
easy axis of the pillar. Using parameters consistent with 
bulk iron, we find that several interesting features emerge 
as the magnetization reverses from an unstable configu- 
ration, saturated opposite the applied field, to the true, 
stable state. Most noticeable is that this highly mis- 
aligned system possesses two distinct decay paths that 
evolve from the same initial conditions. Figure ^ is an 
example snapshot of a slowly reversing decay path, at a 
time when the average z component of the system 
magnetization corresponds to a relatively fiat region of 
the free-energy landscape. As seen in other simulations 
with applied fields along the easy axis of the pillar )S the 
end caps display a region of high vorticity due to pole 
avoidance at the end faces of the pillar. This is shown 
in the left part of Fig. ^ With an applied field nearly 
parallel to the easy axis, the flux lines continue down the 



volume of the pillar, mostly parallel to the easy axis. The 
highly misaligned case discussed in this paper, however, 
exhibits magnetic flux lines that penetrate faces parallel 
to the yz plane near the mid-section, as shown in the 
right part of Fig. 

Many nanomagnets have lateral dimensions on the or- 
der of the exchange length and can be sufficiently rep- 
resented by the Stoner-Wolhfarth mechanism of coher- 
ent rotationi Larger systems that span several exchange 
lengths such as the one discussed in this paper, however, 
are not adequately described by coherent rotation and 
consequently require the application of full micromagnet- 
ics simulations to allow the proper dynamics to emerge. 
To investigate the dynamics of these nanopillars under 
the infiuence of an applied field, we use a micromagnetic 
simulation performing numerical solutions of the Landau- 
Lifshitz- Gilbert (LLG) equation. •^•'^ For our model, the 
pillar's real dimension of 10 nm x 10 nm x 150 nm are 
mapped onto a 6 x 6 x 90 computational lattice, and the 
LLG equation is evaluated at each site in the presence of 
a local field, H{fi). The net saturation magnetization 
Ms is introduced in the LLG equation. 



dt 



1 + a 



H{rl) 



a 



M{r-l) X H{n 



(1) 

as a constant magnitude of every unit magnetization vec- 
tor M{t\) in the lattice. Also included in the equa- 
tion are the constant electron gyromagnetic ratio 70 — 
1.67 X 10^ Hz/Oe and a phenomenological damping pa- 
rameter, whose value we take to be a = 0.1. The local 





FIG. 1: Streamline rendering of a nanopillar before revers- 
ing to the stable state. Mz is denoted by the color of the 
streamlines: +2 is dark gray, —2 is light gray (red/green for 
color online). The applied field represented by the arrow in 
this right-handed coordinate system is oriented 75° from the 
— z-axis, parallel to the zx plane (pointing to the left and out 
of the page). The left part shows the top of a pillar; the right 
part, the same pillar further down the easy axis. 



field iiifi) that determines the dynamics at each site i is 
a smxi of all participating fields, including applied, dipole, 
exchange, and a stochastic thermal field with a Gaussian 
distribution that obeys the fluctuation-dissipation theo- 
rem. In this simulation the applied field was set to 3160 
Oe at an angle of 75° to the easy axis, which is near the 
coercive field for the slow decay mode. The temperature 
was set to 20 K. 

Initially, the nanopillar was allowed to relax in the 
presence of an unreversed field before sinusoidally revers- 
ing the field 180° to drive the reversal process. Measure- 
ments were taken at regular intervals during the simu- 
lation as the LLG equation was temporally integrated 
using an Euler integration scheme. 



II. RESULTS AND CONCLUSIONS 

The first evidence for the existence of two separate de- 
cay modes in the nanopillar comes from the distribution 
of the lifetimes r, determined by the time it took the 
z-component of the total magnetization, Af^, to reach a 
value of —0.75. Lifetimes of decay processes such as this 
one are strongly determined by the shape of the free- 
energy landscape that the system evolves through. For 
free-energy barriers that are much larger than the tem- 
perature ksT), an exponential distribution of the 
lifetimes should be observed. Regions with negligible 




FIG. 2: Cumulative distribution of the lifetimes in nanosec- 
onds. The separation in a fast and a slow mode is clearly 
seen. 
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FIG. 3: Plot of the M^My phase plane for the slow (gray) 
and fast (black) mode. Simulations begin close to center of 
the graph, evolve through the metastable region, and end near 
the top right corner. 



barriers lead to a Gaussian distribution of the lifetimes, 
corresponding to perturbations around a deterministic 
behavior. In Fig. [3 the cumulative distribution for all 
our 17 simulations is shown. The graph is clearly di- 
vided into two regions. The first region corresponds to 
a rapid increase in the distribution for times less than 
about 2 nanoseconds, and the second to a much slower 
increase that is spread across tens of nanoseconds. For 
the fast switches, evidence presented below suggests there 
is no free energy barrier. In that case, we expect the fast 
switching times to have an approximately Gaussian dis- 
tribution. It is not clear what form the distribution for 
the slower switches has. 

Another useful method that reveals the different re- 
versal modes consists of recording the total magnetiza- 
tion (which we normalize to be unity when all spins are 
aligned) of the nanopillar at regular time steps and con- 
structing a phase plot of the evolution. Two examples. 
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FIG. 4: Growth and shrinkage probabihties distinguishing 
the two decay modes in the simulation. The location and 
depth of the metastable well in the free energy depend on the 
reversal mode. Upper panel: fast mode. Lower panel: slow 
mode. 

with T = 1.6 and 17.6 nanoseconds, respectively, are 
shown in Fig.|21for the x and y components of the total 
magnetization. Both configurations begin their evolution 
near the middle of the plot and evolve through their re- 
spective transition states, ending up close to the top right 
corner. There is a difference in the regions where these 
two examples spend most of the time before reversal, and 
this difference is seen for all simulations. This behavior is 
consistent with a local minimum in the free energy for the 
slow mode, and a much shallower (or nonexistent) min- 
imum for the fast mode. However, from the projective 
dynamics evidence to be presented next, it is not clear 
that a true minimum exists, even for the slow mode. 

Projective dynamics^"^'^ gives information contrasting 
the two modes in this nanopillar simulation. By selecting 
a slow variable, growth and shrinkage probabilities are 
accumulated in bins along that variable and can be used 
to further elucidate properties of the free-energy land- 
scape. Specifically, the difference between the growth 
and shrinkage probabilities is a measure of the shape of 
the free energy along that variable. In this instance, Mz 
was chosen for binning, and transition probabilities were 
constructed based on the number of times Mz jumped 
to an adjacent bin. The results of the projective dynam- 



ics analysis are shown in Fig. 0] For the slower reversal 
mode, the growth and shrinkage probabilities practically 
coincide in the region around Mz « 0.7. This indicates 
that the free energy is almost fiat in this region and that 
there is no strongly preferred direction of motion for Mz ■ 
In comparison, the growth probabilities always exceed 
the shrinkage probabilities for the fast mode, indicating 
a tendency for the magnetization to continuously evolve 
toward the stable state. Since neither of these graphs 
exhibits an area where the shrinkage probability exceeds 
the growth probability, it is possible that a significant 
barrier in the free energy does not exist. It is also pos- 
sible that Mz is not a good choice for the slow variable, 
e.g. it may not coincide well with the reversal path. The 
statistics that were collected to construct the projective 
dynamics plots were sorted based on the lifetime of the 
particular run. Due to the approximate classification of 
the modes, it is likely that there is some mixing between 
the data in these plots, although not enough to change 
the general features observed here. 

In conclusion, finite-temperature simulations were per- 
formed on iron nanopillars in a strongly misaligned ap- 
plied field. Sophisticated visualization tools were used 
to give a detailed picture of the magnetization configura- 
tions during the decay. Two distinct reversal modes were 
found, despite identical initial conditions. Three differ- 
ent methods were used to characterize the differences be- 
tween these two modes: the distribution of the lifetimes, 
phase plots of the magnetization transverse to the easy 
axis, and projective dynamics. Still, our results do not 
yet provide a very clear picture of the shape of the free- 
energy landscape. Fully characterizing the free energy 
will be important for designing nanomagnets that switch 
in a reliable manner under given experimental conditions, 
and will be a topic of future research. 
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